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ables, this method saves over an order of magnitude in computer time over previous methods. 
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1. Introduction 

Overlap fermions [jl|] are unquestionably the theoretically cleanest formulation of lattice QCD, 
because they alone exactly fulfil the Ginsparg-Wilson relation [§] and satisfy an exact chiral sym- 
metry on the lattice However, they are also unquestionably the most expensive formulation. 

In [Q], I explored the possibility that the lattice overlap Dirac operator is connected to the 
continuum by a block renormalisation group transformation; this was motivated by consideration 
of their obedience to the renormalisation group derived Ginsparg-Wilson relation. The construc- 
tion contained within that paper was incomplete, because it did not address how to construct the 
Yang-Mills action in a similar way. I will finish this work in a future paper by using a coordinate 
transformation of the gauge fields, and I briefly outline the method in section |23] of this work. Once 
the suitable transformation is constructed, and the correct form of the Yang-Mills action found, then 
it will be possible to transfer from one lattice spacing to another, including the continuum, through 
the renormalisation group. 

This opens up the possibility of running lattice simulations at two length scales: a fine lattice 
for the gauge fields, and a coarse lattice for the fermions. The connection between the coarse lattice 
and an equivalent action on the fine lattice will be provided by renormalisation. The advantages of 
this approach are obvious: if we use a scale factor of s between the fine and coarse lattices, then we 
save at least a factor of s 4 in computer time for the expensive parts of the calculation. However, the 
physics would be identical on the coarse lattice and the fine lattice, which itself can be extracted 
if necessary at an additional cost, but will only be required for certain observables which require 
a finer resolution (for example in studies using heavy quarks). The challenge is linking the coarse 
action with a fine action. One possibility is re-weighting, but the re-weighting coefficients will be 
of order e a 6" b v+ - (the leading order a A V term can be absorbed into the Yang-Mills action), which 
are likely to be unstable at larger volumes without fine tuning of and the terms hidden by the 
ellipsis. A better approach is to find some method of interpolating between the coarse action and 
the fine action; and here I will present one possible way of doing this. 

These proceedings are intended only as a very preliminary description of the proposed method. 
A full proof of the method's validity, numerical results, and a discussion of whether this is indeed 
beneficial will follow in a subsequent paper. In section ||[ I outline the underlying renormalisation 
group theory. In section ||, I outline the proposed Monte-Carlo routine; In section f|, I describe the 
numerical implementation which I will use in my tests, and, in section BL I conclude. 

2. Renormalisation group and overlap fermions 
2.1 Gauge Action 

The multi-grid HMC algorithm presented in section ||| requires that the action is written entirely 
in terms of the Dirac operator (I shall justify this statement in the follow-up paper). This includes 
the Yang-Mills term. I note that in the continuum, if D = y^(e^ — A^), then 

-ia^F^ = (y 5 D) 2 - ^Tr 5 ( 75 D) 2 , (2.1) 
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where the trace is only over the spinor indices. Therefore, the Yang-Mills action can be written as 

S g = |TrF M 2 v = ^Tr[( 75 D) 4 ] - i T r[( 75 D) 2 Tr s [( 75 D) 2 ]]. (2.2) 

This can obviously and easily be extended to any lattice formulation of the Dirac operator, although 
for any Dirac operator which is not ultra local, including the overlap, the numerical implementation 
requires a stochastic estimate, either following the method proposed in [g], or by using a Gaussian 
heat-bath and a rational approximation of the exponential. The idea of using the trace of a function 
of the Dirac operator to simulate the gauge action was first proposed in 

2.2 Renormalisation of fermion fields 

A block renormalisation group transformation of the fermion fields can be defined in terms of 
blocking matrices B 1 and B 1 (which must be functions of the gauge fields) and a (which must 
be independent of the gauge fields) 

Zo= fdU { dxif^dW^e-^-^^f^x 

7=0 J J 

J rf V ^^ e -J^" 1 (^ -^' , CS (,> )- 1 )«(^ -( BW )- 1 ^ ) oc 

Zi = n fdU I d^ ) d^ ) e- s '^-^^^ )+J (2.3) 

7=0 J J 

For simplicity, I will subsequently suppress the flavour indices i and j. For a = °°, we may write 

BD Q B =D { J =Tr[logZ>i - logD ] 

xjf^xjfoB' 1 Yi= B ~Wo- (2-4) 

A possible choice of blocking is 

B 1 =(D0 1/2 D 1/2 ; fi- 1 =D Y '\D^I\ (2.5) 

which leads, if D\ is the lattice overlap operator and Do the continuum Dirac operator, to the 
Ginsparg-Wilson relations, 

7l D 1 +D i7r = 0, (2.6) 

with 

yl,r = y 5 { ±y 5 (»i-»i)( 1_ 4 J ' (2 - 7) 

and an associated ^ ^-symmetric lattice chiral symmetry. 

In [Q, I introduced a parameter £ to control the 'latticeness' of the fermion fields. The Dirac 
operator D\ was re- written in terms of the exponentials, £ 4 e~£E*i Fm-^mI, where « is a lattice site and 
x a position in the continuum, so that as £ — > oo these exponentials become Dirac 5-f unctions and 
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a lattice theory is reproduced. While it is difficult to visualise the theory at the target of £ = °°, at 
finite £ everything remains well-defined. I showed that, for overlap fermions, the blocking remains 
valid as the limit £ — > oo is taken. A full discussion and proof can be found in [Qj. 

In this work, I consider a similar blocking from a fine lattice overlap operator, D , to a coarse 
lattice overlap operator, D\, each with an underlying fine gauge field Uq and E/j respectively. I write 
the blocking, the Dirac operator, gauge fields and fermion fields in terms of a parameter , which I 
shall use to interpolate from the coarse lattice to the fine lattice. We may write 

D^=B^[U^]D [U^]B[U^] (2.8) 

At E, = 0, B = 1 and D remains the original fine Dirac operator; at E, = 1 is equivalent to a 
coarse Dirac operator. The fermion field transformation may be written as 

A /„ ? Xlogto./]^' — — -i/ £ °4rlog[D f ,]^' m 

^ =g 2J0 3|7 61 5 ]i ^. y =^0 e ^ 5 (2.9) 

and an infinitesimal blocking is 

W^ = n + ^^og[D^Mn- (2-10) 

Thus the Jacobian for this infinitesimal change in the fermion fields together with the change in the 
gauge action is 

logdet/- S[S g ] = TrS<^ [S^D]) , (2.11) 
where S| [D] = log[Dg] - S g [Z^] 

2.3 Gauge field flows 

We may define a differential with respect to the gauge fields as 



itT a jj 



For a gauge field transformation given by, 

j^J^=M^[U^]U^, (2.13) 

the Jacobian, Jjj, of the infinitesimal transformation satisfies [0] 

Trlog/fy = iS«^yn^ iM [E/^]). (2.14) 
Thus, the solution to the equation 

Tr[^(n° M )+^(5)n^+a^] =q, (2.15) 

where d? excludes the dependence on E, within the gauge fields, and Q is independent of the gauge 
fields, completes the renormalisation group link between Do and D\ (and equally the renormali- 
sation group construction of the overlap action). It can be shown that a solution for exists for 
all £ as long as S and d^S are analytic, and integrating over this solution will give the required 
transformation of the gauge fields. Finding this solution is rather more problematic. However, it 
is possible to map from the coarse gauge field to the fine gauge field stochastically, and for an 
importance sampling method, that is good enough. 
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3. Multi-grid HMC 

The detailed balance equation reads, 

p[u8^uZ]w c [o,ut] =p[u£ ^ug\w c [o,u§] =p[u^ut]p[ul^u?]p[u?^u8]w c [o,uZ], 

(3.1) 

where, in t/J, the subscript refers to £ and the superscript to Monte-Carlo time and 

N f 

W c [^U^ = Y\dct[D i Me- s ^. (3.2) 

!=1 

If each of the three steps satisfies detailed balance, then it is clear that the complete update of the 
gauge field satisfies detailed balance. For the transfer from the fine to the coarse lattice, I introduce 
the momentum field IT according to the Gaussian distribution g~( n+A ?) , where 

A^ = (d^(S))-\d^. (3.3) 

This requires that (d" fl (S)) is invertible (I shall leave the proof until the later work), and that the 
Dirac operator is constructed so that is independent of the gauge and pseudo-fermion fields at 
t; = and t; = 1. We then proceed by integrating along the equations of motion using a reversible, 
area conserving procedure, such as the Omelyn integrator, along a trajectory T, 

j^U =i(U+A^)U; =d^S[D^] - (3.4) 

which conserves E = (Yl^ + A% ) 2 + TrS[D^ ] . IT now is a conjugate field and does not depend on the 
gauge field. The Jacobian from the gauge field update cancels the Jacobian from the momentum 
update. Reversibility can be maintained using an iterative procedure. I estimate the trace using 
pseudo-fermions generated according to a heatbath and a rational approximation R[D^\ ~ e s ^ D ^ 
obtained using the Remez algorithm, so that / d<j>d<j>^e~^ R ^ = detR 1 = e ^ TlS ^ D \ The detailed 
balance equation for the transformation from [Ui] to [Uq] reads (if the integration is exact) 

P[U <-Ui]Wc[Ui] = J dllrdUiJ d^d^e-^+ A ^\-^ R ^8{[U Q ,U ( {\-T[U i ,U i \) 

= J dnfJUMd^e-^+^e-^Pd+SipuTli] -T- l [U ,Uo]) 

=P[U l ^U ]Wc[Uo}. (3.5) 

We can correct for the error from the numerical integration either by re-weighting or an extrapo- 
lation in the integration step-size. Similarly, the integration from Uo to U\ can be shown to satisfy 
detailed balance, although for practical purposes this step will not be required, because each coarse 
configuration can be generated from the same thermalised fine configuration. 

We run enough coarse trajectories to gain an independent configuration, and, if required, gen- 
erate a fine configuration using the procedure above. For many quantities, the gauge field asso- 
ciated with the fine Dirac operator will not be needed: calculations can proceed with the coarse 
Dirac operator and non-perturbative renormalisation will take care of the rest. Only if the finer 
lattice resolution is needed for a particular observable will it be necessary to perform the costly 
conversion to the fine lattice and the equally costly calculation of the observable. 
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4. Numerical implementation 

The initial kernel operator, Kq, is defined as a Wilson operator on the fine gauge field: 

(*oW = 75 ^ -£[(!- 7;.)^/^ W + (1 + r,.)5,-^(/;^ - A)] j • (4-1) 

Using a scale factor of three, I define a blocking operator, B l K for the kernel in terms of the parameter 
% and additional parameters C l t (which can be tuned to minimize the force and ensure stability when 
constructing the fine gauge), 



(B^yy [q = 1 1 u z>jJ [(i - - (i - ef^n, 

n 6 (I - \yp - n P l) 6 (I - \y'p - n n l) . c 4 - 2 ) 

where G 1,4,7, .. ., £ x „ a path of links between x and n (constrained within the hypercube), and 
the corresponding path ordered product of links. For a coarse kernel operator ^ = (B^KoB^ 
the Fourier transform (in the free theory) takes the form iy^ sin^ + (1 — cos(p^)) + liy ll sin(3p jU + 
(1 — cos(3p /J )) + m for some i. For i > 1 this introduces fermion doublers. I have chosen to avoid 
these doublers by including a second term with different coefficients, so that 

Ks = (B° K ?K B° K - a(B l K ^KoB l K + P(B K )^B° K - ^{B l K )^B l K - y(l - (1 - £ 2 ) 2 )(1 " 8 x , n 8„y), 

(4.3) 

where the parameters a, j8, j3' and y need to be tuned to avoid doublers. A badly tuned set of 
parameters is usually obvious from the eigenvalue spectrum of the kernel and overlap operators. 
Bk is constructed so that 9^Bk = at t; =0 and t; = 1, as required by the algorithm described in 
the previous section. The overlap operator at mass jj, is D[E,\ = (1 + (1 — /^/ssign^). 

Construction of the HMC algorithm for the coarse Dirac operator is trivial. I am currently 
testing it on a small PC-Cluster, using a 4 3 8 coarse lattice to generate 12 3 24 gauge field ensembles, 
and comparing it against a 12 2 4 trajectory generated using normal methods; timing the coarse 
multi-grid HMC over trajectories of length 1 and the non-multi-grid comparison over trajectories of 
length 0.02 and extrapolating to a length 1 trajectory. The multi-grid code gave, on average, a factor 
of 107 gain over the original routine. There was, however, an overhead of a factor of three from a 
normal 4 3 8 run, which was mostly caused because both the overlap and Kernel operators were very 
poorly conditioned compared to the small lattice counterpart, and partly from the additional work 
needed to differentiate the gauge fields. Obviously, the cost of constructing the fine lattice has yet 
to be included in this estimate. It is to be hoped that with proper tuning of the kernel operator 
that this cost will be no more than that of a single length 1 trajectory on the larger lattice, meaning 
that overall the gain will be roughly proportional to the autocorrelation length. 

5. Discussion 

I have discussed whether it is possible, using overlap fermions, to generate fine gauge fields 
ensembles using a coarse Ginsparg-Wilson Dirac operator. The method will not be valid for other 
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lattice Dirac operators because the underlying renormalisation group theory breaks down. Small 
scale tests show that this method gives a two orders of magnitude gain over my previous HMC 
algorithm, although with additional factors not yet included in that cost; the most important of 
which is the cost of extracting physics from the gauge field ensembles. It may reasonably be asked 
what use is there in generating the ensembles quickly using this method if one then has to apply an 
overlap operator on an exceptionally large lattice to extract physics. For most quantities it will also 
be possible to do the physics on the coarse lattice: the coarse and fine lattice observables will be 
linked to each other by renormalisation, and the renormalisation constants can be calculated non- 
perturbatively as usual. The number of observables which require the fine lattice will, hopefully, be 
quite small. It should also be possible to adapt this method to allow an anisotropic formulation on 
the fine lattice. Why then use this method at all, since it gives a considerable overhead compared 
to the usual approach of a coarse Dirac operator on coarse gauge field? Because the locality of the 
overlap operator is improved (which I have confirmed numerically), it will (if correctly constructed) 
have smaller lattice artefacts, and it will have a greater sensitivity to the topology of the fine gauge 
field, and, of course, the fine gauge field and larger lattice is available for those observables which 
require it. 

However, the question of how much of a gain, if any, is achievable by this approach is not really 
answerable until it has been fully tested. I hope to present numerical results in the subsequent paper. 



Acknowlegments 

I am grateful to the support of the Deutsche Forschungsgemeinschaft grant DFG FOR-465, 
and for many useful discussions with Andreas Schafer. I am also grateful to the hospitality of the 
Kivali Institute of Theoretical Physics in Beijing during the workshop "Lattice Quantum Chro- 
modynamics" in July 2009. Numerical calculations were performed on the PC Cluster Juli at the 
Jiilich Supercomputer Center in Germany. 

References 

[1] R. Narayanan, H. Neu berger, Chiral fermions on the lattice, Phys. Rev. Lett. 71 (1993) 325f-3254. 
arXiv:hep-lat/930801 1; R. Narayanan, H. Neuberger, Chiral determinant as an overlap of two vacua, 



Nucl. Phys. B412 (1994) 574-606. arXiv:hep-lat/9307006; H. Neuberger, A practical implementation 



of the overlap-Dirac operator, Phys. Rev. Lett. 81 (1998) 4060^062. arXiv:hep-lat/9806025 



[2] P. H. Ginsparg, K. G. Wilson, A remnant of chiral symetry on the lattice, Phys. Rev. D25 (1982) 2649; 
[3] M. Liischer, Exact chiral symmetry on the lattice and the Ginsparg-Wilson relation, Phys. Lett. B428 



(1998) 342-345. | arXiv:hep-lat/9802011 



[4] N. Cundy, A renormalisation group derivation of the overlap formulation, Nuclear Physics B824 



(2010)42. arXiv:0903.5521 



[5] B. Joo, I. Horvath, K. F. Liu, The Kentucky Noisy Monte Carlo algorithm for Wilson dynamical 



fermions, Phys. Rev. D67 (2003) 074505. arXiv:hep-lat/01 12033 



[6] I. Horvath, A framework for systematic study of QCD vacuum structure. II: Coherent lattice QCD 
arXiv:hep-lat/060703 1 ; 



[7] M. Luscher, Trivializing maps, the Wilson flow and the HMC algorithm arXiv:0907.5491; 



7 



